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Global attractivity and positive almost periodic 
solution of a multispecies discrete mutualism 

system with time delays 
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Abstract — In this paper, we consider an almost periodic 
multispecies discrete Lotka-Vol terra mutualism system with 
time delays. We first obtain the permanence and global 
attractivity of the system. By means of an almost periodic 
functional hull theory and constructing a suitable Lyapunov 
function, sufficient conditions are obtained for the existence of a 
unique strictly positive almost periodic solution which is 
globally attractive. An example together with numerical 
simulation indicates the feasibility of the main results. 

Index Terms — Almost periodic solution, Mutualism system, 
Discrete, Global attractivity 


I. Introduction 


Recently, investigating the almost periodic solutions of 
discrete and continuous population dynamics model with time 
delays has more extensively practical application value(see 
[1-19] and the references cited therein). In this paper, we are 
concerned with the following multispecies discrete 
Lotka-Volterra mutualism system with time delays 


Xi{k+ 1) = Xi(k) exp di(k) - bi(k)xi(k - a*) + 
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Xj(k - Tfj) 


<*«(£) + x i( k - Tij) 


, i = 1 , 2 , •• • ,n , ( 1 . 1 ) 


where (a/k)}, {/?/£)}, {Cij(k)} and {d^(k)} are bounded 
nonnegative almost periodic sequences such that 


0 < a- < a.i(k) < a“, 0 < b\ < bi(k) < 6 “, 

0 < c‘j < Cij(k) < cfj, 0 < d[j < dij(k) < d“ , 


z,j = 1,2, , n,j ^ i, k E Z. For any bounded sequence {f(k)} 

defined on Z, f u = sup/(A:),/ / = inf f(k). 

k £ z 

By the biological meaning, we will focus our discussion on 
the positive solutions of system (1 . 1). So it is assumed that the 
initial conditions of system (1.1) are the form: 

Xi(0) = > 0, v?i(0) > 0, 

0 G N [— r, 0] = { — r, — r + 1, 0}, 

t = max {<7j, Tjj}. (1.2) 

To the best of our knowledge, this is the first paper to 
investigate the global stability of positive almost periodic 
solution of multispecies discrete Lotka-Volterra mutualism 
system with time delays. The aim of this paper is to obtain 
sufficient conditions for the existence of a unique globally 
attractive almost periodic solution of the systems (1.1) with 


initial condition (1.2), by utilizing an almost periodic 
functional hull theory and constructing a suitable Lyapunov 
functional and applying the analysis technique of papers [3, 
12, 13]. 

The remaining part of this paper is organized as follows: In 
Section 2, we will introduce some definitions and several 
useful lemmas. In the next section, we establish the 
permanence of system (1.1). Sufficient conditions for the 
global attractivity of system (1.1) are showed in Section 4. 
Then, in Section 5, we establish sufficient conditions to 
ensure the existence of a unique strictly positive almost 
periodic solution, which is globally attractive. The main result 
is illustrated by an example with a numerical simulation in the 
last section. 


II. Preliminaries 

First, we give the definitions of the terminologies involved. 
Definition 2.1 ([20]) A sequence x: Z^R is called an 
almost periodic sequence if the e-translation set of x 

E{s,x} = {t £ Z : | x(n + r) — x(/z) |< e.Vn G Z } 

is a relatively dense set in Z for all e > 0; that is, for any given 
e>0, there exists an integer /(s) > 0 such that each interval of 
length /(e) contains an integer x ^ E{e, x} with 

| x(n + r) — x(n) |< s, Vn E Z. 

x is called an e-translation number of x(n). 

Definition 2.2([21]) Let D be an open subset of Rm, f : Z 
xD — »R m . f(n, x) is said to be almost periodic in n uniformly 
for x £ D if for any e > 0 and any compact set S c D, there 
exists a positive integer /=/(e, S) such that any interval of 
length / contains an integer x for which 

1/(7! +T,x) — f(n y x) | < £, V(7!,t) E Z X S. 

x is called an e-translation number of f(n,x). 

Definition 2.3( [22]) The hull of f, denoted by H(f), is 
defined by 

H(f) = {g(n,x) : lim f(n + t*,, x) = g(n,x) uniformly on Z x 5}, 

k—+oo 

for some sequence {x k }, where S is any compact set in D. 

Definition 2.4 Suppose that X(k) = (x^k), x 2 (k), • • • , 
x n (k)) is any solution of system (1.1). X(k) is said to be a 
strictly positive solution in Z if for k^Z and i = 1, 2, • • • , n 

0 < inf Xi(k) < sup Xi(k) < oc. 

kGZ 

Now, we state several lemmas which will be useful in 
proving our main result. 

Lemma 2.1( [23]) (x(n)} is an almost periodic sequence if 

* 

and only if for any integer sequence { k t }, there exists a 
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* 

subsequence { k { } c { k i } such that the sequence {x(n+ k t )} 

converges uniformly for all n £ Z as i — mx>. Furthermore, the 
limit sequence is also an almost periodic sequence. 

Lemma 2.2( [24]) Assume that sequence {x(n)} satisfies 
x(n) > 0 and 

x(n + 1) < x(n) exp{a(n) — b(n)x(n)} 

for n^N, where a(n) and b(n) are non-negative sequences 
bounded above and below by positive constants. Then 

lim supx(n) < yyexp{o ri — 1}. 

n— H-oo O 

Lemma 2.3( [24]) Assume that sequence {x(n)} satisfies 

x(n + 1) > x(n) exp{a(n) — 6(n)j , (n)}, n > Aq. 


lim sup x(n) < x* 

n— M-oo 


and x(N 0 )>0, where a(n) and b(n) are non-negative sequences 
bounded above and below by positive constants and N 0 ^ N. 
Then 


lim inf x(n) > min 

n— ► -foe v 




Xi{k - (7 i) < Xi(k)ex p { - (Ti[a[ - 6“ (Mi + £■)]}. (3.7) 

Substituting (3.7) into system (1.1), for k>K+^, it follows 
that 

Ti(k + 1) > Xi (k ) exp |flj - 6- exp { - - 6- (Mi + s)]}xi(k) J. (3.8) 

Thus, as a direct corollary of Lemma 2.3, according to (3.4) 
and (3.8), one has 

lim inf Xi(k) > min{mj ie ,mj 2£ }, (3.9) 

k—t+oo 

where 


m 


m 


He = -j£ exp {ai[a\ - + e)]}, 

= rn il£ exp |a| - b“ exp { - + e)]}Mi |. 


Letting e— >0, it follows that 

lim inf Xi(k) > — min{mij ,raj 2 } = zn* > 0, 

k—t+oo 2 

where 

m h = j ^ ex P - Vi M i)}, 

m i<2 = m ix exp ja- - 6“ exp { - a l (a\ 


(3.10) 


- b?Mi)}Mi 


III. Permanence 


In this section, we establish the permanence result for 
system (1.1). 

Theorem 3.1 System (1.1) with initial condition (1.2) is 
permanent, that is, there exist positive constants mi and M\(i = 
1, 2, • • • , n) which are independent of the solutions of system 
(1.1), such that for any positive solution (x^k), x 2 (k), • • • , 
x n (k)) of system (1.1), one has: 

nii < lim inf Xi(k) < lim sup Xi(k) < Mi , i = 1,2,*-' , n. 

k-*+oo fc-H-oc 

Proof. Let (x^k), x 2 (k), • • • , x n (k)) be any positive solution 
of system (1.1) with initial condition (1.2). From the first 
equation of system (1.1) it follows that 

n 

Xi(k + 1) < Xi(k) exp {a.i(k) + Y2 Cij(k)} 

n 

< Xi(k) exp {a“ + c ij }• ( 3 -l) 

By using (3.1), one could easily obtain that 


x i {k-a i )>x i {k)exv{-(j i (a“+ Y c“)}. (3.2) 

Substituting (3.2) into the first equation of system (1.1), it 
follows that 


( n n y 

<+ Y c «“ - b \ ex p { - a *( a “ + Y 


(3.3) 


Thus, as a direct corollary of Lemma 2.2, according to (3.3), 
one has 


lim sup Xi(k) < -j exp{(a" + Y + 1) - 1} - M*. (3.4) 

fc->-+oo t i 

For any small positive constant 8 > 0, from (3.4) it follows 
that there exists a positive constants K > 0 such that for all k > 
K and i = 1, 2, • • • , n, 


Xi(k)<Mj+s. (3.5) 

For k > K+cr„ from (3.5) and system (1.1), we have 

Xi(k + 1 ) > Zi(fc)exp { di(k ) - bi(k)xi(k - Vi)} 

> Xi(k) exp {a - - b x (M l + £)}. (3.6) 


Thus, by using (3.6) we obtain 


Then, (3.4) and (3.10) show that system (1.1) is permanent. 
The proof is completed. 


IV. Global attractivity 


In this section, by constructing a non-negative Lyapunov- 
like functional, we will obtain sufficient conditions for global 
attractivity of positive solutions of system (1.1) with initial 
condition (1.2). We first introduce a definition and prove a 
theorem which will be useful to obtain our main result. 

Definition 4.1 A solution (x^k), x 2 (k), • • • , x n (k)) of 
system (1.1) with initial condition (1.2) is said to be globally 
attractive if for any other solution (x*x(k), x* 2 (k), • • • , x* n (k)) 
of system (1.1) with initial condition (1.2), we have 
lim (x*(k) — Xi{k)) = 0, i = 1,2,* • • ,ra. 

k — >-+oo 


Lemma 4.1 For any two positive solutions (xx(k), 
x 2 (k), • • • , x n (k)) and (x*i(k), x* 2 (k), • • • , x* n (k)) of system 
(1.1) with initial condition (1.2), we have for k > 2x 


In 


if(t+l) *?(i) 

71 

+ Y. «y(*Ki(*)- 

J= 1j¥« 


--6t(*)[r,(fc) -<(*)] 

X](k-T tj )-r*(k-T tJ ) 
[dij(k) + xj{k - Tij)][dij{k) + Tj ( k 



k-l , n 

-h : \k) Y, |h,( 5 )-^(*)]^s )[%(«)— Y 

s=k—o, ' 




Xj(s-Tij)- Xj(s-Tij) 


+ i 1 (s)B i (s) 2 Cijisjdtjis)- 

KuF) + W - T '\J )]^j(s) + X*{s - T ( j)] 

- bi(s)[ii{s -a,)- x*(s - ff,-)] 


(4.1) 


where 

Ai(s) = exp j 0,(.«) [<!,(.<?) 


xUs — Tij) ,1 

bMx ‘ (s “ ffi) + }' 


B,(s) = exp 


| <Pi («)[«.(») ~ f>i(8)Xi(8 -<t ( ) + 
+ (1 - <fii(s))[a,(s) - b,{s)x*(s 


Y c (s) - T 'J ) 1 

Jb A ) d ii {s) + x J (M-r ij y 

\ V" i \ r j( s ~ T tj) 



( 4 . 2 ) 


€ (0. 1), i = 1,2.- • • , n. 

Proof. For i = 1, 2, • • • , n, we can have from system (1.1) 
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ln _ ln = ln _ ln <%+_!) 

*i(* + l) x *{k ) M fc ) 

Xj(k — Tij) 


= ai(k) - bi(k)xi(k - (Ti) + 53 


3=1,3'/* 


d y (A:) +Xj(fc - nj) 


x*(A'-r 0 ) 


a, W -M*)x*(*-.,) + E ^W <;yW + I .( fc - Tu )j 

3=1,3/* J 3 J 


_ 53 c ij(^) 

3=1,3'/* 


Xj(A’ - Tij) 


X*(k-Tjj) 


_dij {k) + Xj(k- Tij ) dij(fc) + x* ( A- - Tij ) 


6i(fc)[x t -(fc - <Ti) - x*(fc - <Ti)] 


+ %(» - r y )]Kj(*) + 4(* - T„)] '"WM*) 

+ bi(k){[xi(k) - Xi(k [x* (k) - x*(k - o-j)]}, 

that is 

+ 1) _ x<(fc) i ^ _ (l .|j , ,. x J'jCA* — Ty) — Xj (A- — r,j) 

x;(k + l) x*(k) 2- 1 [d,j{k) + Xj (k - r^ld^k) + xj(k - r tJ )] 

-bi{k)[Xi(k) - rr* (A*)] + bi(k){[ii(k) - ii{k - <r,)] - [x*(0 - x,*(A* - a,)]}. 

Since 

[Xi(k) - Xi(k - (Ti)\ - [x*(k) - X* (k - <x<)] 
k - 1 k-1 

= Y [Xi{s + 1) - Xi(s)\ - Y K(s+t)-2:-(s)] 


(4.3) 


S=k — (Ti 

k-1 


s=k—Oi 


= Y {[^i(s + t) + t)] - [Xi(s) -x*(s)]}, 

S=k — <Ti 

and for k>2x 

[Xi{s + 1) - X* (s + 1)] - [Xi(s) - X* (s)] 

= x i (s)exp[a i (s)-6i(s)xi(s-<j i )+ Y c *' 3 '( s )' 


(4.4) 


3=1,3'/* 


dij(s) +Xj(s-T lj ) 


] 


Xjis-Tij) 


x*(s)exp[a i (s)-b i {s)x*(s-<Ti)+ Y Ci ^ dr(s \ + x *i s _ Tr) 

3=1, 3V* *3V / j\ *3/ 


] - [Xi(s) -x*(s)] 


[*«(*) “®i(*)]| exp[a i (s)\-b i (s)xl(s-CT i )+ Y c *i( s ) 




d i:j (s) + x*(s 


3=1, j/i ,JV ' ' ^ 


r]- 1 } 


+ x i (s){exp[a i (s)-6 i (s)x,(s- ( T i )+ ^ c y (s) ^) + Xj (s - r 0 - ) - 


3 = 1,3'/* 
n 


exp[a i ( S )-6 i ( S ) I *( S -^)+ £ 

3 = 1, jV* J 3' J 


Using the Mean Value Theorem, we get 

Ms + 1 )-x*(s+ 1)] - [ Xi{s ) - X* (s)] 


= [Xi(s) - x*(s)]Ai(s) 


X*{s-Tij) 


a.(s) - 6 ,(.*(. - a,) + JE «M Ms) + i;(s _ Ty)J 


+ X , 


( S )B,( S ) [dy („) + 4(« - ry)][dj(>) + 4(» - r«)l 


-6 f (s)[Xi(s - ffi) - x* (s - <Ti)] 


(4.5) 


here A,-(s), Z?,(s) are defined by (4.2). Then from (4.3)-(4.5), 
we can easily obtain (4.1). The proof is completed. 

Theorem 4.1 Assume that in system (1.1) with initial 
condition (1.2), there exist positive constants /?,(/ =1, 2, • • • , n) 
and tj> 0 such that 


PiEj - Y Pj F ij > *?> * = 1, 2, • • • , n, 

3=1,3#* 

where 

2 


(4.6) 


B, = min{Sj, jjr - »f> ~ - <*<»?*?(«? + *f«< + E ^jr*). 

* 3=1, 3V* 


•^*3 = "^"(1 + a jMjbjBj). 
*3 


(4.7) 


Then for any two positive solutions (x^k), x 2 (k), • • • , x n (k)) 
and (x*!(k), x* 2 (k), • • • , x* n (k)) of system (1.1) with initial 
condition (1.2), we have 

lim (x*(k) - Xi(k)) = 0, i = l,2,---,n. 

fc-^+oo 

Proof. Firstly, let 

Vn(k) = | \nxi(k) — In a;* (A:) | . 

From (4.1), we have that for k > 2x, 


3=1 , J 


k- 1 


+b,(k) y {i*i(*)-*r(*)i^i(*)[«*w+^wi*r(*-»i)i+ Y 

*=k-Oi 3=1 

+ |x,(s)|iii(s) Y l®3 1*3 ) Xj(s Tij ) | + 6j(s) |x, (s <7| ) X^ (s <7; ) |1 \ . (4.8) 

j=l . j¥i a 'A*) JJ 

Since 

Ti(fc) - X*(k) = e lnx *W _ e ln *JW - k(k) \n( Xi (k)/x*(k)), i = 1,2,- • • ,n, 

where £,(k) lies between Xi(k) and x*i(k), i = 1, 2, • • • , n, it 
follows that 

| In (xi(k)/x*(k)) - bi(k)[xi{k ) - xj(fc)]| 

= |ln(M*)/*?(Ar)) -& t (%(A-)ln(x,(fc)/x*(fc))| 

= I ln( Xi (k)/xt(k))\ - (~ - |^y - 6,(fc)|) |® 4 (*) - x* (*)|. (4.9) 

By Theorem 3.1, there are constants Mi> 0, and a positive 
integer k 0 such that for k > k 0 , 0 < Xi(k), x*i (k) < M„ i = 1, 
2, • • • , n. Then from (4.8) and (4.9) we can obtain that for k > 
k 0 + 2x , 

AV " W 5 -(m- - Id*) - b ‘ {k) \h (t) - «*>i + i i§, - T “ ] - ^ - r «)i 

3= ljV* 1 

+b,(k) x; |^o)[«,(s) + a/a(s)+ y 


3-1 


+ MiS.(s) V ^f(|x3 ( a -r 0 )-x;(*- r.j)! 

+ A/iB i («) & *(s)|x.(s - <Ti) - x*(s - o-j)| j. 

Secondly, let 

v,,a-)= e e adi±M| Ij( .)- i; (.)| 


(4.10) 


3=1,3^* B=k — T, 

k- 1+ff, Jb — 1 


53 53 {‘4t(«)[o,(«) + 4/,6,(t‘) + 53 , 1^ x,(u) - x;(?i)| 

s—k u=*-a, k 3=1 / 


+ A/,Z?,(u) 53 ^ 7 ( ' ») 1^4^ ~ r «3 )~ Jj(«-T 0 )| 


3=1 jV* 


+ A#iS<(u)6t(«)|£i(u ~ ^i) ~ ^*(« ~ 

By a simple calculation, we can obtain 

AV', 2 (i')= E E 


(4.11) 


3=1, 3V« 


d, j (A - + r,j) J j=1 ^ | . d,3 ( k ) 


+ (.4, <*)[*(*) + A/ i 6 i (*)+ E -Vj|4^]hi(*-)-i?(t)| 

^ 3=1, 3^i ^ ’ 


+MiBi(k) 53 

3=1 jV« j 

4 l'+ ff i 

+jw*)m*)M*' 0 *) -«*-*)!} E 6 '< s > 


fc— l 


s=A*+l 

n 


“In (k) Y {A(«»)[fli(i«) + Mi6.(w)+ 53 4/ ? 2 ^|]|x,(u)-x*(u)| 

u=k-a t 1 3=1jV* ,jW 

-:U,B,(u) 53 T^yl z i( u - T y) ~ “ T *f)| 

3=1,3/* 


+ A/ ( Bj(u)6,(u)|ii(u -<7j) - x*{u -<7j)j|. 

Thirdly, let 


(4.12) 


A’-l l+Tii+O-j 


Viz(k) = Mi EEL Wi(* + ru) ^|l ^ - »;(0| 

i-utu-k-T,, s=i+Tij+i a w + 

A:— 1 /+2CTj 

+ Mi E E 6 i (s)B,(/ + <J i )b,(; + a, )|ij(/) -I*(/)|. 

l — k — Oi s = i + <7j + l 

Then we can derive 

n I, \ A+Tjj+<7. 

AV i3 (k) = M i E Biik+n /f Mft)-3#)| E 6 '( s ) 

3=1,3#* *=fc+T«+l 

fl /i \ 

-MiBi(k) Y ^7n\ x j( k - T ij)- x *j(k-rij)\ Y b ^ s ) 

3=1,3#* s=fc+l 

j 

+ M i B j (fc + a j )Ml + <*i)hW-irW| E 6 '( s ) 

S=fc+(X i -fl 
k+(Ti 

-M i Bi(fc)6 i (fc)|x j (i:-<T0-x*(r--<Ti)| E 

s=/c+l 

Now we set 
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Vi(k) = Vn(k) + V i2 (k) + V i3 (k), i = 1, 2, • • • ,n. 

Then from (4.8)-(4.13), we have that for k > k 0 +2x , 


AVi(k) < - 


1 


1 


P c ij{ k + x ij) 


6W 

ft If / 7 \ 

+ Ai(k)[ai(k) + Mibi(k) + Y ‘-&M.] Y 6 i ( s )| Ii (fc)-x*(fc)| 

j=1j¥* y ' 




S — fc-f 1 


Cjj (fc T Tjj) 




+Mt y w+Tn) fy—y y ms)|i#)-i*(*)| 

lJ lJ S=fc+Tji + 1 


A'+ 2 cr , 


+ + Vi)bi(k + Gi) Y “ x i ( k ) |- 

Now we define a Lyapunov-like discrete functional V ( k ) by 


n 


v(k) = ym(k). 


i=l 


It is easy to see that V(k 0 +2x)<+oo. Calculating the difference 
of V(k ) along the solution of system (1.1) with initial 
condition (1.2), we have that for k > k 0 +2x, 

n f k+ 20i 

i=l V. t\ / i s=Ar+<Tj+l 


-.«*) [«,(*) + A/AW + Y PPlT 


k-cr, 


] Y w 


s=fc+l 


- E ffrfrf- E wm® E 

j-lj* W + ’W j.L.jA Oj'l* + 7 J.I„ t + r , i + | J 


'j 




k-T^+aj 


<£ (ft ["**{*!. U W - - Mwfc+W + E ^)] 

»=i l 1 j=1j¥« a ' J 


J=U*t ij 1 

n n 

= ~Y ft®*- E WW*)-#) 

t=l 

n 


<-,Eh«-*;«|, 


«=i 


where E l and /q,- are defined by (4.7). 

Then we have that 

k k n 

Y [v{p+i)~ v(p)] < -t} y Y Mp) - x i ^ I ■ 

p=k 0 +2r p=k 0 +2r i=l 

which implies 

k n 

V(k+l) + T) Y X!Mp) _x i(p)| <U(fc 0 + 2r). 

p=k 0 +2r i= 1 

That is 

k " V(k 0 + 2r ) 


H 5Z |^.(P) -<(P)| < 


and then 


p=k 0 +2r i= 1 


+oo n 




Y ] < 


At — /lqT^t i — 1 


V{ko + 2 r) 

n 


< +OG, 


71 


which means that ^ lim Y |^i(^) — £*(A.’)| = that is 
lim (x t(k) — x*(k)) = 0, i = 1, 2, • • • , n. 

fc— H-oo 

It means that (x^k), x 2 (k), • • • , x n (k)) is globally attractive. 
This completes the proof of Theorem 4.1. 


V. Almost periodic solution 

In this section, we will study the existence of a globally 
attractive almost periodic sequence solution of system (1.1) 
with initial condition (1.2) by means of an almost periodic 
functional hull theory and constructing a suitable Lyapunov 
function, and obtain the sufficient conditions. 

Let { S m } be any integer valued sequence such that S m —>00 
as m— >oo. According to Lemma 2.1, taking a subsequence if 
necessary, we have 

Gi(k + 6 m ) — > a*(k),bi(k + £ m ) — » b*(k),Cij(k + <5 m ) — > c*j(k),d.jj(k + S m ) — > 


d*j(k ), i,j = 1, 2, • • • , n, j ± i, as m — > 00 for k G Z. 

Then we get a hull equation of system (1.1) as follows: 


Xi(k + 1) = Xi(k) exp < a*(k) - b*(k)xi(k - er,) + 


c * ( m x j( k ~ T >j) 1 

i= t^ tjU d* j (k) + X j (k-T ij )f' 


i= 1,2,--- ,n. (5.1) 


By the almost periodic theory, we can conclude that if 
system (1.1) satisfies (4.6), then the hull equation (5.1) of 
system (1.1) also satisfies (4.6). 

By Theorem 3.4 in [26], we can easily obtain the lemma as 
follows. 

Lemma 5.1 If each hull equation of system (1.1) has a 
unique strictly positive solution, then the almost periodic 
difference system (1.1) has a unique strictly positive almost 
periodic solution. 

Theorem 5.1 If the almost periodic difference system (1.1) 
satisfies (4.6), then the almost periodic difference system (1.1) 
admits a unique strictly positive almost periodic solution, 
which is globally attractive. 

Proof. By Lemma 5. 1 , we only need to prove that each hull 
equation of system (1.1) has a unique globally attractive 
almost periodic sequence solution; hence we firstly prove that 
each hull equation of system (1.1) has at least one strictly 
positive solution (the existence), and then we prove that each 
hull equation of system (1.1) has a unique strictly positive 
solution (the uniqueness). 

Now we prove the existence of a strictly positive solution 
of any hull equation (5.1). By the almost periodicity of {a*i 
(k)}, |b*i(k)}, |c*ij(k)} and |d*ij(k)}, there exists an integer 
valued sequence {x m } with x m — ► 00 as m — ► 00 such that a*i (k + 
7m) -► a*i (k), b*. (k + T m ) ^ b*. (k), c*,j(k + i m ) -> c-^k), d-./k 
+ Tm) — ► d*ij(k), as m— >oo for k ^ Z. Suppose that X(k) = (x^k), 
x 2 (k), • • • , x n (k)) is any solution of hull equation (5.1). By the 
proof of Lemma 2.2 and 2.3, we have 


nii < liminf x i (fc) < lim sup Xi(k) < i — 1, 2, • • • , n. (5-2) 

fc-H-oo fc-^+oo 

And also 


0 < inf Xi(k) < sup Xi(k) < oo, i = 1,2, ••• , n. 
fcez+ fcGZ + 

Let 8 be an arbitrary small positive number. It from (5.2) 
that there exists a positive integer k 0 such that m-s^x^k)^ Mi 
+ s, k > ko, i = 1 , 2, • • • , n. Write X m (k) = X(k + x m ) = (x im (k), 
x 2m (k), • • • , x nm (k)), for all k > k 0 + x - x m ,m £ Z + . We claim 
that there exists a sequence {yi(k)}, and a subsequence of {x k }, 
we still denote by {x k } such that x im (k) yi(k), uniformly in k 
on any finite subset B of Z as m— >oo, where B = {ai, a 2 , . . . , 
a p }, a h ^Z(h = 1, 2, . . . , p) and p is a finite number. 

In fact, for any finite subset BcZ, when m is large enough, 
x m + ah — x>k 0 , h = 1, 2, . . . , p. So 

m t - e < Xi(k + r m ) < M{ +e, i = 1 , 2 , • , n, 


that is, (xi(k+x m )} are uniformly bounded for large enough m. 

Now, for ai ^B, we can choose a subsequence {x (1) m } of 
{x m } such that {xi(ai+x (1) m )} uniformly converges on Z + for m 
large enough. 

Similarly, for a 2 ^ B, we can choose a subsequence {x (2) m } 
of {x (1) m } such that {Xi(a 2 + x (2) m )} uniformly converges on Z + 
for m large enough. 

Repeating this procedure, for a p ^ B, we can choose a 
subsequence {x (p) m } of {x (p_1) m } such that {xi(a p +x (p) m )} 
uniformly converges on Z + for m large enough. 
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Now pick the sequence {T (p) m } which is a subsequence of 
{x m }, we still denote it as {x m }, then for all k^B, we have 
Xi(k+x m )— »yi(k) uniformly in k^B, as m— >oo. 

By the arbitrary of B, the conclusion is valid. 

Combined with 

im{k + 1) = Tim{k) exp s Gj (k -f- T m ) bj (k -j- T v k)Xi rn (k g t) -b 


lim \x*(k) - y*(h)\ = 0, i = 1,2, 

k->— oo 


(5.5) 


n 


Define Q = ^t^iQi , where 


1=1 


1 ^ TijCV , r . " MjCV 

' J l J , ' 2 lu f 4 u /„u , ?i r lu , \ J l J 


Gi = -+ E 

in. - £ — ' 






+ a?f>j‘[.l?(a? + M,»r + X 


f"- 


, )+«,£?( X # + &?)] 




n 


c£(fc + T m ) 


%jm{k Tij ) 


+ AW?B?( E 


U T- ■<* 
'V c y 


~h Tm) ~h Xjm (k Tij ) J 

* = 1, 2, • • • , n, 




+ (Tib?), i = 1,2, — ,n. 


gives 


y,(A- + 1) = yi(k)ex p <j a*(fc) - h](k)y,(k - a t ) + 

Vjik-Tij) 


E <#(*) 


<*<#(*) + iy(*-fy) 


, i = 1,2, • • • ,ra. 


We can easily see that Y(k) = (yi(k), y 2 (k), • • • , y n (k)) is a 
solution of hull equation (5.1) and m, - e < y,(k)< M t + s, i = 1, 
2, • • • , n, for k ^ Z. Since 8 is an arbitrary small positive 
number, it follows that m, < y*(k) < M„ i = 1, 2, • • • , n, for k ^ 
Z, that is 

0 < inf yi (k) < sup yi(k) < oo, * = 1,2,--- , n. 
k ^ Z z 

Hence each hull equation of almost periodic difference 
system (1.1) has at least one strictly positive solution. 

Now we prove the uniqueness of the strictly positive 
solution of each hull equation (5.1). Suppose that the hull 


Let 8 be an arbitrary small positive number. It follows from 
(5.5) that there exists a positive integer K > 0 such that 

\x*(k)-y*(k)\ < ^,k<-K,i = 1,2,... , n. 

Therefore, for k<-K, i = 1, 2, • • • , n 

v ai k ) < E tf 1 m <k _ y *^ P ^ 


%J 




+ A? (a? + Mib? + E '—JT 1 ) ™ I x i(p) - y*{p ) I 

y p - 


n cY 


+ MiBf E ~jT \ x j&) ~ VjiP ) I + MiBfb? max |a£(p) - y,*(p)| 


j=l,j^i 13 

X I #+^y[ J 4!-« + M i 6v + £ ! • A/.»r! X c f 

3=1 Jjti lj 3=1,3^* y y 

\2 /lu\ 2 nu 


VS(fc) < E 


71 T' -/? — 
1 l J c l J 


J 1 UCU 1 

<M,a,f>“B“( £ 

j=i jv* y ^ 


max |x*(p) - y*(p)| + (6“) max |x*(p) - y*(p)| 

t><L- J J p <k 


equation (5.1) has two arbitrary strictly positive solutions It follows from (5.3) and above inequalities that 

V*(k) <^(3iQi^=£, k<-K , 


(x*i(k), x* 2 (k), • • • , x* n (k)) and (y*i(k), y* 2 (k),- • • ,y* n (k)). 
Like in the proof of Theorem 4.1, we construct a Lyapunov 
functional 

n 

V(fc) = X^(VuW + ^(*) + V' , 3 (fc)), fce Z, (5.3) 

i=l 

where 


1=1 


n 


fc-1 


Tij ) 

j=l,j^j s=k-T lj ^ lJ 


vsw = E E 


dij( s + Tij) 




A: — 1-j-o'i 


fc-1 


n 


+ E ^ A( w )[ a t( w ) + Mibi(u) + E 


s=k 


u—s—a 
n 




AljCjjju) 

dij(u) 




+ M,Bi(u) Yi f4^Tl*J(“- r «)-#J(“- T «)| 




+ A/i^i(w)6i(w)|x*(w - (Tj ) - y*(iz - 

^SW = ^i E E E W-Wa + T-yJ^^llarJW-yJWI 

j=l,J7^i l—k—r ij S=l+T ij+ 1 + 7 '3> 


k — 1 i+2<Tj 


l = k — (Ti S=i+<Tj + l 

Calculating the difference of V*(k) along the solution of the 
hull equation (5.1), like in the discussion of (4.14), one has 

n 

A V(k) < -r,Y k*(*) - w?(*)l. k e z. 


so limV*(fc)=0 • Note that V*(k) is a non-increasing 

k — ^ go 

function on Z, and then V*(k)=0. That is x*i (k) =y*i (k), i = 1, 
2, • • • , n, for all k^Z, Therefore, each hull equation of 
system (1.1) has a unique strictly positive solution. 

In view of the above discussion, any hull equation of 
system (1 . 1) has a unique strictly positive solution. By Lemma 
2. 2-2. 3 and Theorem 4.1, the almost periodic difference 
system (1.1) has a unique strictly positive almost periodic 
solution which is globally attractive. The proof is completed. 

VI An EXAMPLE AND NUMERICAL SIMULATION 

In this section, we give the following example to check the 
feasibility of our result. 

Example Consider the following almost periodic discrete 
Lotka-Volterra mutualism model with delays: 

ii(k + 1) = xi(A')cxp / 0.025 + 0.005 sin(v / 2A*) - (1.0075 - 0.0025 cae(\/3fc))xj (A: - 1) 


+(0.03 - 0.005 cos(\/2fc)) 


x 2 (k- 2) 


1.02 + 0.005 sin( v/2 k) + j •■,(k — 2) 


(5.4) 


(0.02 + 0.000 sin(yAik))- 


x 3 (k - 3) 


i= 1 


From (5.4), we can see that V*(k) is a non-increasing function 
on Z. Summing both sides of the above inequalities from k to 
0, we have 


4.03 + 0.005 cas(y/3k) + x 3 (* - 3) }' 
r 2 (k + 1 ) = x 2 {k) exp 1 0.035 + 0.005 cos( xftk) - ( 1 .0025 + 0.0015 sin(v / 2/fc))x 2 (ifc - 2) 


+(0.02 - 0.00 lsi»( v/2A*))r— 


ii(k - 1) 


0 n 


vYYi - »?(«)! ^ - V '(k + !)> * < 0. 

q=k t=l 

Note that V*(k) is bounded. Hence we have 


5.03 + 0.(X)G eos( \/3k) + x, (k - 1 ) 

+(0.025 + 0.005 cos(\/5fc)) J:<(A ~ 2) 1. 

5.01+0.01 sin(v/2A ) + j* 3 (A* - 2) J 


( 01 ) 


0 n 


q= — o o i=l 


which implies that 
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x 3 (k + 1) = i 3 (A)exp jo.02G - 0.006 sin( y/&k) - (1.0035 + 0.0025 cos(>/2*))*3(* - 1) 

+(0.03 + 0.005 cos(v4tt*)) Xl(A '~ 2> 

4.08 - 0.004 sin( v/3A-) + x x (k - 2) 

+(0.02 - 0.006 sin(v^Ar)) ~ 3> 1. 

4.06 + 0.008 cos(y/2k) + x 2 (k - 3) J 

By simple computation, we derive 
Mi « 0.0261, M 2 « 0.0392, M 3 « 0.0283, 

E\ % 0.0347, F 12 ~ 0 . 0177 , F 33 « 0.0152, 

E-2 ~ 0.0358, F-21 ~ 0.0136, F 23 ~ 0.0168, 

E 3 « 0.0407, F 31 « 0.0145, F 32 « 0.0183. 

Then 

E\ — F 12 — F 13 ~ 0.0018 > 0.001, 

E 2 - F 21 - F 23 « 0.0024 > 0.001, 

F 3 — F 31 — F 32 ~ 0.0021 > 0 . 001 . 

Also it is easy to see that the condition (4.6) is verified. 
Therefore, system (6.1) has a unique strictly positive almost 
periodic solution which is globally attractive. Our numerical 
simulations support our results (see Figs. 1-3). 




FIGURE 1: Dynamic behavior of X)(k) of system (6.1) with 
the initial conditions (x^k), x 2 (k), x 3 (k)) = (0.015, 0.019, 
0.011) and (0.034, 0.022, 0.028), k= 1,2, 3, 4 fork e [1, 100] 
and k £ [500, 550], respectively. 




FIGURE2: Dynamic behavior of x 2 (k) of system (6.1) with 
the initial conditions (x^k), x 2 (k), x 3 (k)) = (0.015, 0.019, 
0.011) and (0.034, 0.022, 0.028), k= 1,2, 3, 4 fork e [1, 100] 
and k G [500, 550], respectively. 




FIGURE3: Dynamic behavior of x 3 (k) of system (6.1) with 
the initial conditions (x^k), x 2 (k), x 3 (k)) = (0.015, 0.019, 
0.011) and (0.034, 0.022, 0.028), k= 1,2, 3, 4 fork G [1, 100] 
and k G [500, 550], respectively. 
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